K-essence scalar field as dynamical dark energy 
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We study an early dark energy (EDE) model as a K-essence scalar field in the framework of FLRW 
universe using an effective parametrization of the state equation as a function of the redshift z with 
the tracker condition during radiation domination, but also demanding an accelerated expansion of 
the universe at late times emulating cosmological constant. We found all the dynamical variables 
of the EDE system. We use the luminosity distances of the SNIA to get the best estimations for 
the free parameters of the model and also, we constrain the model using primordial abundances of 
light nuclei in BBN theory. We summarize the necessary conditions to achieve BBN predictions 
and the accelerated expansion of the universe at late times. 
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I. INTRODUCTION 

Recent observations of the luminosity distances of the SNIa |19| reveal that the expansion of the universe is ac- 
celerated and there is an unknown matter-energy contribution as about 70% of the critical density, which is smooth 
and has negative pressure. In order to explain this phenomenon, it has been proposed many plausible solutions: the 
cosmological constant A which is related to the vacuum energy of the quantum fields [T], Quintessence fields (with 
the state equation ui = ^=constant), Kessence, Taquionic fields, frustrated topological defects, extra-dimensions, 
massive (or massless) fermionic fields, galileons, effective parametrizations of the state equation, primordial magnetic 
fields, holographic models, etc. All these proposals have been used to model this contribution predicted by the 
Friedmann equations in the framework of the General Relativity. Also, there are other possible options as Modified 
Gravity, where the accelerated expansion effect is geometric and not as a matter-energy form. 

Nowadays, the current paradigm is the ACDM model, which is so far the best fitting to the present observations, 
even though the conceptual problems that persist with the nature of A. 

On the other hand, primordial abundances of the light nuclei that were formed during Big Bang Nucleosynthesis 
(BBN) are well observed and quantified with astrophysical methods, specially the mass fraction of the *He. Actually, 
there is an extended theory of BBN, introduced by Alpher, Bethe and Gamov [2ii3j, and it has been developed many 
numerical codes that resolve the Boltzmann equation for each isotope including the cosmological background (based 
on Kawano code). There is a good agreement between the predicted abundances and the observations, however, the 
main problem remains with Li. Despite the efforts that have been made, including corrections on the cross sections 
of this element, the discrepancies are not negligible. 

With the aim to enhance the primordial abundances calculated from BBN and also, give a plausible explanation of 
the accelerated expansion of the universe, we have proposed a model of dark energy which has a non-null contribution 
at early times to increase the Hubble radius during radiation domination era and influence the Boltzmann equations 
that determine the evolution of the light abundances. All the conditions that allow us to describe the early dark 
energy are achieved with a K-essence scalar field, which is characterized by its state equation that overcomes the 
attractors defined by a dynamical system. 

In the first section, we expose the main conditions imposed to the scalar field and we resolve the dynamical system 
that appears from the cosmological assumptions. At the second part, we propose the effective parametrization of the 
state equation, determinate the best estimations of the free parameters of the model using the luminosity distances 
of the SNIa from the Union dataset and derive analytical expressions for the dynamical variables of the K-essence 
system. In the IV section, we test our model with some standard proves and summarize the main differences between 
this model and ACDM. 

Taking into account that there is a non-negligible energy density of the K-essence field during radiation era, we 
compute the BBN abundances, including the scalar field degrees of freedom in Hubble parameter, following the 
derivation and discussion made by Bernstein |25J and using the available codes for Nucleosynthesis. Finally, we find 
the values oi Qb and rjB predicted by our model as a result of the BBN calculations. 



II. K-ESSENCE SCALAR FIELD AND THE DYNAMICAL SYSTEM 

The K-essence scalar fields appeared at the late 90s with the K-inflation model proposed by Armendariz-Picon 
[71|S]. However, the idea was extended to describe a dynamical dark energy contribution, taking into account that 
this field can track during radiation domination epoch and also, it could avoid the fine-tuning of the initial values 
of the field and its velocity. These features are well known for the system, so the challenge is resolve the evolution 
equations of the field in the FLRW spatially flat universe and without cosmological constant. 

The lagrangian of the K-essence scalar field is given by: 

p(X,0)=X(<^)i(X), (1) 

where X = — |V"0Vq,(/) is the kinetic energy of the field and v its velocity v = ^ — \/ —IX > 0. Rewriting ([T]) in 
terms of v, the lagrangian has the following form : 
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p{v,<l>) = K{cf>)Q{v). 

The most general action that describes the K-essence field in a cosmological plasm is given by: 



(2) 



R 

2k^ 



(3) 



with — SttG and 5b is the action of the background matter. The signature that is used is {— 1,+1.+1,+1} The 
equations of motion of the field come from the variation of the lagragian with respect to the field (j): 



1 dp d f Idp 

V dv dv \v dv 

From ([3]), it is obtained the energy-momentum tensor 



^^(■jH) I V - — 
dv d(j)dv d4> 



(4) 



df^(j)d„(l)-p{(j),v)gf^t,. 



(5) 



Since the scalar field can be described as a perfect fiuid, the energy density and the pressure are defined by: 



= i^(0) [v^-Q], 



(6) 



p^ = K{^)Q{v). 

In addition, the adiabatic velocity of sound for the K-essence field is given by: 



(7) 



vQ" 



(8) 



where / denotes derivative with respect to the velocity of the field v. gives relevant information of the stability of 
the perturbations associated with the K-essence field. 

In order to resolve the ^ and find an explicit form of ^ and ([t]), there have been suggested many alternatives: 
fixing a specific function of (f>{t) or v{t) [H], making redefinition of the field to face a modified lagrangian Q{v) [23\, 
considering pure kinetic K-essence model |10fll31 flS] or imposing Slow roll conditions on K[<j)). 

However, we want to resolve the complete dynamical system defined by the Friedmann and the continuity equations 
for non-interactive fiuids in the cosmological background: a matter (or radiation component) and the K-essence scalar 
field. 



H 



(9) 



- = -7r((l + 3w0)/90 + (1 + 3w„i)/9,„), 
a 6 



(10) 



(11) 



Reexpressing the equiations (|9|, ( 10 1, (111 in terms of the velocity and the dimensionless variable F = pf-j"p (matter 
energy density fraction): 



dv 



{lnK)_^v f Kp^ 



1/2 



(12) 
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dF 3k 



(13) 



In addition, it must be fulfilled the following conditions: 



dp4, _ (1 + uj^) _ 



(14) 



dv ^ V \ CI 
It is assumed an asymptotic behaviour for the function i^((/)) [TS] : 



(15) 



X(0) = ^i§^, lim/fo(0) = O. 



(16) 



To reach an attractor during radiation domination epoch, the dynamical system {w(0), i?((^)} must fulfill the ansatz: 



(17) 



where A{(j)), B{(j)) — > monotonically for — )■ cx) (or equivalently, v — Vrad and F ^ 0). The ansatz implies the 
following physical conditions on the system: 



U!,n = W0(wo), P0(0) 7^ 0, C^(wo) > '^n 



(18) 



^'M = 1 - 



1 + w,, 



^0. 



(19) 



There is also a De-Sitter attractor given by the condition R ^ 0, which guarantees the existence of a accelerated 
expansion at late times: 



The last ansatz entails a condition on the state equation on the vicinity of the De-Sitter attractor: 



(20) 



1 - UJ^jVs) 
1 + UJ^{Vs) 



>0, \uj^{Vs)\<l. 



(21) 



III. EFFECTIVE PARAMETRIZATION OF C04, 



In order to get a general solution for the dynamical system defined by (14 1, (151, (18 1, (19) and (211, it has been 
proposed an effective parametrization of the state equation from z < 10^^ given by: 



uj^{z) = 



4/3 



/ \ ■''1. 



1, 



(22) 



where m is factor that modules the transitions between the attractors, z^ is a redshift in matter domination epoch 

z ~\~ z 

defined by Zfi — — and z^, the redshift where the De-Sitter domination -accelerated expansion- begins. 
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The parametrization (22) respects all the conditions previously mentioned, hence it is possible to resolve the 



functions related to the K-essence lagrangian. 
The energy density of the field ^ results: 



integrating ( 23 1 , it is obtained: 



(23) 



p = po • (1 + zY 



)? 



4/m 



Po ■ f{z). 



with 



/(a) = exp 



_3 , ^ (i±M^,,; 



(^) 



i/r. 



Meanwhile, the fraction of the dark energy density £7^ = 1 — F 



(24) 

(25) 
(26) 



f^</.o • /(a) 



f^0o • /(a) + i^mo • a 



(27) 



The formal solution of ( 22 1 will be obtained with the best estimation of the free parameters of the model {il^o , m, } . 
For this reason, we use the luminosity distances of the SNIa from the survey Supernova Cosmology Project with z > 0.8 
to minimize the function 



[pi - p{Zi)] 



(28) 



with the distance modulus given by the expression p = m — M = 5{logiodL{z) — 1) and the luminosity distance of 
our model: 



dUz) 



c(l + z) r dz' 



Ho Jo B{zy 



(29) 



B{z) = {n^J{z'-m, z,) + (1 - + z'ff'^ 



(28 1 must be resolved together with the constraints: 



E^^^^ I Hi - p{zi) \ / dp{z.i;^^^,m,z^) 

i—1 ^ ^ 



= 0, 



(30) 



where x — Q^g^m, z^,. In addition, we have imposed other 2 conditions for the parameters: the deceleration param- 
eter has to be zero at z.^,, therefore q{z.^,) — 0, but also f{z ^ zbbn ~ 10^) has to overcome the maximum value 
at the Primordial Nucleosynthesis to contribute with some relativistic degrees of freedom and enhance the predicted 
primordial abundances. 



Resolving (28 1 simultaneously as (30 1, it is found the following values for the free parameters: 
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Figure 2. Function f{z): dark energy density evolution. 
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In the figure 1 is plotted the evolution of the state equation uj as function of a. The field emulates radiation during 
this epoch and then it evolves to the next attractor: De-Sitter. Here is clear the tracker behaviour imposed in the 
dynamical system. In addition, it is shown that the field is relaxing to the asymptotic ACDM model for late times 
(during De-Sitter attractor.) 

Figure 2 shows the behaviour of the the function f{z) that characterizes the dark energy density evolution in the 
model. During radiation domination epoch, the field scales as radiation p oc (1 + z)** until Zgg. After that, p has a 
complex behaviour which guarantees the second attractor will be reached. At this point, the K-essence scalar field 
evolves in the De-Sitter attractor and its state equation goes asymptotically to —1 (as cosmological constant). 

The figure 3 displays the luminosity distance for the model compared with the predicted by ACDM. The shift 
between the curves shows that the luminosity distance is upper than the associated with the ACDM model, be- 
cause the matter density today predicted by our model is higher than the second one (compared with WMAP-7 
{17^0, ^mo} = {0.734 ± 0.029, 0.266 ± 0.029} [? ]. 

However, for low redshifts the luminosity distance grows linearly independent on the model at this regime. 

On the other hand, figure 4 displays the evolution of the dark energy density fraction for the model with the 
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Figure 3. Luminosity distance as a function of z. 
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Figure 4. Evolution of the dark energy density fraction as a function of z. 



estimations of the free parameters (III) 



Up to now, we are interested to resolve the dynamical variables of the K-essence according to the effective 
parametrization we have proposed. The first one of these quantities is the adiabatic velocity of the sound C^, 
which has a behaviour defined by equation (31 1. 



m 



Zd 



1 + z 



(W0 + 1). 



(31) 



Using (22 1, we complete the evolution of C^. The behaviour of the adiabatic velocity is plotted in the figure 5. 



fulfills the condition ( [l8| in the radiation domination epoch attractor. Whatsoever, there is not a straightforward 
physical interpretation of the adiabatic velocity during De-Sitter attractor, because this value implies a complex Young 
module of the plasma perturbations. 

It might be that the field is unstable when < 0, however, it is necessary to compute the evolution of the field 
perturbations and the time of the stability condition to conclude something in this respect. 
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Figure 5. Adiabatic sound velocity evolution within the model. 



On the other hand, the velocity of the field (the time evolution of the field) has a solution in terms on the scale 
factor given by: 



Wo • (1 + z) 



l+Zd 
1+z 



1 



4/m— 1 



(1 + ^d)" + 1 



(32) 



where vq is the field velocity today. 

Finally, the non-canonical term of the action results: 



(/ \m , 4/m+l 

(1 + =,) +1) 

That is the predicted behaviour for the non-canonical term of the lagrangian in terms on z. More interesting is that 
Q{v) evolves as radiation in the first attractor and then acquires a more complicated dependence during the transition 
between the attractors, ensuring the continuity of the function. 



IV. STANDARD COSMOLOGICAL PROVES 



Age of the universe with this model 



The age of the universe according to our model is given by the expression: 



'"-irJ! (o.3i»-3 + a69/(«))V^ '^° = ^ 

This is just an approximation of the age of the universe; however, this is an excellent result taking into account that 
the parametrization was figured out for z < 10^^. 



On the other hand, the existence of this scalar field implies that the universe evolves faster than in the standard 
model, because for a more negative w, more accelerated is the expansion and older is the universe for a given Hq. 
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Matter inhomogeneities evolution 



When the Meszaros equation is resolved within the model with the K-essence scalar field, one solution for the 
inhomogeneous modes grow with the Hubble factor de Hubble. Those are the modes which maintain their amplitude 
after cross the horizon. 

There are many effects on the density of perturbations: it is suppressed the linear growth co{z) < —1/2, with respect 
ACDM model, where the growth is proportional to the scale factor a; this suppression raises with higher uj at matter 
domination epoch and entails an earlier beginning to the dark energy domination epoch (accelerated expansion) . 

Otherwise, considering the CMB anisotropics spectrum, it is possible to verify that for values of to > 1.0, there 
is formation of the first two acoustics peaks, then the model has a strong influence during radiation and matter 
domination epochs, which not correspond with the observations. On the other hand, for to values lower than 0.5, the 
spectrum goes faster to ACMB, and they are rejected because do not scale as radiation in that domination but tend 
W0 ~ —1 for a « a*. 



The acoustic peaks correspond to the modes which at the CMB decoupling epoch were in the maximum compression 
(odds peaks) or rarefaction (even) and their position is really sensitive with the state equation of the dark energy. 
Actually, the first peak depends on fl^^^ monotonically and taking into account that our model predicts a lower value 
for this parameter compared with the ACDM, hence the first peak undergoes a shifting to lower multipole moments. 
However, there is a degeneracy that must be broken down in this phenomenon: the existence of this kind of model 
implies a rise in the number of baryons fi^, that also entails a shifting to lower multipole moments of the peaks. 
Furthermore, the larger fig is, the higher would be the first peak. 

Therefore, if increases its value within the model, the relative distance between the peaks decreases, because 
the mass associated with the baryons makes that the oscillations occur faster. 



Statefinder parameters 



In order to distinguish between our model and ACDM, |14| has proposed a test using Statefinder parameters, 
defined by: 



r = l + ^fl^u;4l+uj^)~^fl^j^, (34) 

When the parameters are computed, we concluded that the K-essence scalar field model is slightly off ACDM, 
because the prediction of the values f2„iOj and are upper than the obtained with cosmological constant. 
However, as we have argued in this section, the model is relaxed and tends to ACDM for redshifts during De-Sitter 
domination epoch. 



CMB shift parameter R 



This parameter measures the shifting of the acoustic peaks from from BAO and it is defined as the comoving 
distance between the last scattering surface and today: 



R={nraHlY'^ j^^ (36) 

with H{z) depends strongly on the model. 

The measured value for this parameter is R = 1.719 ± 0.019 fI7^, meanwhile the numerical calculations made with 
our model give Real = 1.75; the convergence to this value is insured for z < 200. 
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V. CONSTRAINS WITH BBN ABUNDANCES 

As it has been proposed in the introduction of this paper, the main goal with the effective parametrization of 
the state equation for the K-essence scalar field is to avoid the need to fine tuning the initial conditions, but also 
reproduce the BBN assess enhancing the predicted abundances. The devise is really simple, we include the field energy 
density as relativistic degrees of freedom during radiation epoch, then compute the Hubble factor and consequently 
the Boltzmann equations for the light nuclei. 

Firstly, we set out a fiducial baryon to photon ratio rj = 6.2 x 10^^ [? ] and calculate the capture temperature that 
deuterium began to be formed via nuclear reactions given by: 



Tc,T, = eo /26 = 0.088MeF. (37) 

For this temperature the deuterium is not longer photodissociate f25j, where e/j is the deuterium binding energy. 
Temperatures lower than Tc--y the whole chains related with primordial nucleosynthesis could be run. 

Following the discussion of Bernstein , we assess analytically the value of the mass fraction Yi^^ and then, using 
FastBBN [HI HZ], Public BBN and BBNreactions jil8|, we compute numerically the light abundances. 

In order to introduce the effective degrees of freedom of the K-essence scalar field, we impose the condition mentioned 
in section III to restrict the value of the m-parameter: we want to achieve the largest contribution of the energy density, 
been subdominant with respect to the radiation energy density. The condition can be quantified as: 



P4,\rad = b ■ Prad 0<b<l. (38) 

We introduce this energy density in the Hubble parameter during radiation era and executed the time capture, the 
neutron's fraction at this time and finally, the mass fraction, obtaining the following results: 





6 = 


Model with a contribution of fa = 0.2 




1.13 


1.2379 




182 


169.8 


x„(tc) 


0.123 


0.125 




0.247 


0.249 



The values in the |V] are compared with the calculated values within the standard model (where 6 = 0, i.e. null 
contribution of the K-essence scalar field) . 

The abundances for the other nuclei are not precisely computed to be reported using this analytical method, 
because they depend on the coupled Boltzmann chains. 

However, it is notable that the predicted values for the ^He is in agreed with the observational boundaries, therefore 
our model is an excellent candidate for dynamical dark energy model in FRWL. 

Using FastBBN |3[T7], Public BBN and BBNreactions [H], we compute numerically the light abundances including 
the relativistic degrees of freedom of the field as 6: 
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Code 


b 


D/H X 10-5 


^He/H X 10-5 


Yp 


^Li/H X 10-1° 


Fast BBN 


0.0 




2.335 


1.546 


0.241 


1.268 


Fast BBN 


5.0 X 10- 


3 


2.345 


1.548 


0.241 


1.261 


Fast BBN 


5.0 X 10- 


2 


2.436 


1.570 


0.245 


1.263 


Fast BBN 


0.1 




2.537 


1.594 


0.249 


1.261 


Fast BBN 


0.2 




2.741 


1.639 


0.256 


1.269 


BBN reactions 


0.0 


4.043 


2.363 


0.243 


1.543 


Public Big bang 


0.0 




1.542 


3.000 


0.242 


4.884 


Public Big bang 


5.0 X 10- 


3 


1.547 


3.002 


0.242 


4.907 


Public Big bang 


5.0 X 10- 


2 


1.5845 


3.015 


0.246 


5.111 


Public Big bang 


0.1 




1.629 


3.030 


0.250 


5.345 


Public Big bang 


0.2 




1.717 


3.059 


0.257 


5.825 


our EDE model 


0.2 






0.249 




m 


0.0 


2.75 ±0.24 


0.93 ±0.055 


0.2484 ± 0.0004 


3.82 ±0.66 


m 


0.0 


2.60 ±0.18 


1.04 ±0.04 


0.2479 ± 0.0004 


4.15 ±0.47 



Moreover, it is possible to assess the baryon to photon ratio from the CMB temperature anisotropies spectrum, 
taking into account the relation: 

In the table 4 are shown the values of the parameters related with BBN (il^ and tje) in terms of b, where b = 0.0 
is the associated value to ACDM and b — 0.2. the maximum contribution of the field during radiation domination 
epoch: 



P4' 






r]B X 10-1° 


& = 0.0 


0.044 


0.02218 


6.20 


& = 0.2 


0.053 


0.02692 


6.73 


m 


0.0449 ± 0.002S 


0.02258 ± 0.00057 


6.190 ± 0.145 



According to the results shown in the table|v] we plot the abundances for the light nuclei including ^iJe, for different 
contribution of the field, implemented as effective relativistic degrees of freedom in the code Public BBN. 
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VI. CONCLUSIONS 



In this work, we have made a general description of the dark energy component as a K-essence scalar field that 
evolves during radiation as a tracker (w^lrad = |) and at late times achieves the De-Sitter attractor uj^ -> — 1, 
emulating cosmological constant A. All the description was made in the hot cosmological plasm including matter- 
radiation components and assuming a spacially flat universe. With these assumptions for the field, we have proposed 
an effective parametrization for the state equation (22), that its free parameters estimations were obtained by 
minimizing the function with the distances modulus of the Type la Supernovae. 



We have rewritten the velocity of the field and the non-canonical kinetic term of the K-essence lagrangian in terms 
of z, obtaining the completed behaviour of the field during the thermal history and this let us say that K-essence 
system includes different classes of Quintessence, when it is considered specific cases of the non-canonical term. 
Another advantage with this formulation is that avoids the Fine-tuning of the initials conditions of the field and its 
velocity. 

As it was expected, during radiation and matter domination epochs the field has upper predictions for some observ- 
ables rimOj ^4>o and (but respecting the observational bounds for them) comparing with ACDM model. However, 
at late times z ~ tends asymptotically the standard model, after it has evolved from the radiation to the second 
attractor. All these results are in agreement with the conditions that have been imposed to resolve the dynamical 
system |16| . 

On the other hand, when it is included the non-null contribution of the field during radiation domination era, the 
Hubble factor is affected, but also the time capture and therefore, the primordial light nuclei abundances, because 
there were more neutrons out of the equilibrium to form '^He by two body-reactions. Actually, the whole reactions 
occurred faster, such that the production of the nuclei are more effective and drives in an upper mass fraction. 

The predicted value for the ^He abundance prediction according to our model is inside the observational bounds, 
then, our model is an excellent candidate to a dynamic dark energy with a subdominant contribution during radiation 
and matter epoch. 

Finally, it is remarkable that the results of this paper can be compared with other kind of models, because the field 
degrees of freedom can be treated as effective degrees of some other component (for instance, a Quintessence scalar 
field or even, sterile neutrinos In fact, we have include the field contribution in the numerical codes in these way, 
and it let us to compare degenerations of different kind of models. 
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